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Abstract. In this work, a lower bound for the ground state energy of the Falicov- 
Kimball model for intermediate densities is derived. The explicit derivation is 
important in the proof of the conjecture of segregation of the two kinds of fcrmions 
in the Falicov-Kimball model, for sufficiently large interactions. This bound is 
given by a bulk term, plus a boundary term of the form a 1 (n)|9A|, where A is the 
region devoid of classical particles and n is the density of electrons. A detailed 
proof is presented for n = 1/2, where the coefficient a±(l/2) = 1CP 13 is obtained, 
for the two dimensional case. Although clearly not optimal in terms of order 
of magnitude, this is the largest explicitly calculated coefficient in the range of 
intermediate densities. With suitable modifications the method can also be used 
to obtain a coefficient for all densities. That is the topic of the last section, where 
a sketch of the proof for n < 1/2 is shown. 
PACS numbers: 71.10.Fd, 71.10.Hf, 71.30.+h. 

1. Introduction 

The Falicov-Kimball model [ |FKj ] was introduced to investigate metal-insulator transi- 
tions in mixed valence compounds of rare earth and transition metal oxydes. Later, it was 
again considered to describe order in mixed valence systems and binary alloys. A review 



of exact results for this model can be found in [GM]. 

The model assumes two kinds of fermions in the lattice f2: classical (infinitely massive) 
'ions' with density n c = iV c /|J)| and electrons with density n e = N e /\£l\. For simplicity, 
the particles are assumed to be spinless (without loss of generality, the spin variable can be 
introduced later). The Falicov-Kimball hamiltonian can be written in the second quantized 
form 

H = - y~] txyalay + U y~] n x w(x), (1.1) 

where at and a x are the fermion creation and annihilation operators for the electrons in 
x, and n x = a\a x . The variable w(x) can be either 1 or 0, according to whether the site 
x is occupied by a classical particle or not. We will assume here € "L d 

For a bipartite lattice $7 = A U B, Kennedy and Lieb HKLfl proved that the ground state 
displays crystalline long range order at half filling (n c + n e = 1). This result illustrates 
the relevance of the model in fundamental problems in condensed matter physics, like 
understanding the formation of crystals and molecules. Also, it expected that the better 
understanding of the Falicov-Kimball model will provide new insights to the Hubbard 
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model. And in the context of the Hubbard model, other fundamental question can be 
addressed, like the existence of ferromagnetism in a system in which the spins are itinerant 
(not localized). 



A long standing conjecture for the Falicov-Kimball model [FF| was that, for sufficiently 



strong interactions, the two kinds of particles should segregate away from half-filling. This 



conjecture was proved in [FLU| where it is shown that the total ground state energy is 
bounded above and below by a bulk term, plus a second term which is proportional to the 
boundary of the region A devoid of classical particles. If E^,n is the ground state energy 
for N electrons, 

e(n)|A|+ai(n)|0A| < E KN ^ e(n)|A| + a 2 (n)|dA|, (1.2) 

where e(n) is the energy per site for a density n = iV/|A| of free electrons in the infinite 
lattice Z d . Therefore, given that the bulk term is fixed for all configurations, lowering the 
energy requires minimizing the boundary, which is accomplished by segregating the two 
species of fermions from each other. 

Also in |FLU | explicit coefficients a\{n) are obtained for low densities of electrons 



n ^ \Sd\/(^) d , where A is the domain devoid of classical particles and \Sd\ is the volume 
of the d-dimensional sphere, whereas a2(n) is determined for all densities. The lower 
bound is obtained by considering first the U = oo case. Taking t xy = 1, the hamiltonian 
acting on a function <p(x) 6 L 2 (A) can be written 

[h\(p\(x) = 2d(p(x) — cp(x + e), (1.3) 

e:x+e€A 

where the sum is over the edges of the lattice. The eigenvalue equations are h\<pj = ej<pj, 
for j = 1, . . . , |A|. Their lower bound is derived from the inequality 

E KN -\k\e{n) > j^-- !(e F -e k ) £ -l-\(b k , Vj )\ 2 dk, (1.4) 

where e k = 2d — 2 cos h- L . Also, the concept of the boundary vector 

h(x) = XdA(x)e-* kx J2 ^ (1-5) 

is introduced. Therefore, the problem reduces to showing that the boundary vector has 
a projection in the subspace spanned by the largest eigenvalues. The mathematical re- 
sults are bounds for the sum of the lowest eigenvalues of the Laplace operator. For the 



continuous Laplace operator, one should refer to [LY]. 

This will be the starting point of our study here. Our goal is to obtain an explicit coef- 
ficient for the boundary term for intermediate densities \Sd\/{2ir) d < n < 1 — \Sd\/(2ir) d . 
We are going to obtain results for U = oo, from which the results for finite interaction can 
be derived (see flFLUH ). O ur main result in this limit is: 



Theorem 1.1. For d = 2 and density of electrons n = 1/2, the ground state energy of the 
Falicov-Kimball model is bounded below by 

E A , N - |A|e(l/2) ^ ai(l/2)|5A|, (1.6) 

where Qi(l/2) = 1(T 13 . 
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2. The boundary term for n = 1/2 

2.1. Projection of the boundary vector. The goal is to prove that the boundary 
vector b k has a projection in the subspace spanned by the eigenfunctions {<Pj}, &j > 2d 
(h\<pj = ejtfj). If we can prove that this projection is proportional to the boundary for 
a subset (of non-zero measure) of the region in fe-space limited by the fermi surface of 
n = 1/2 (ef = 2d), the boundary term can be calculated. 
Expanding bk in terms of the eigenfunctions of h A we have 

-£|foA)| 2 (e;-2d)+2 \(<Pj,h)\ 2 (e j -2d) = Y,\(<Pj>h)\ 2 \e j -2d\ 

3 j:e 3 >2d j 

2 (gj - 2d) 2 \\(h A - 2d)b k \\ 2 



2d 

Therefore 

Y. Kw.WI'fe-M) > llf*-JOM' + fa.(*A-aoW t (21) 

j:ej>2d 

and 

£ ite.«i 2 > [K^z|W + (M^M s/w . ( , 2) 

j:ej>2d 

Suppose we can find fc such that e k = 2d and \\(h A — 2d)bk\\ 2 ^ a|<?A|, for some 
constant a. For such k, we have to consider the two possible cases: 

• (b k ,(h A -2d)b k ) > 

• (b k , {h A - 2d)b k ) < 

We should only be concerned with the second case, where the negative contribution 
from the second term could cancel out the boundary term. We claim that for k' = 
k + (ir,ir, ...,7r), e k , = 2d, \\(h A - 2d)b k >\\ = \\(h A - 2d)b k \\ and (b k ,,(h A - 2d)b k ,) = 
-(b k ,{h A -2d)b k ) > 0. 

Indeed, if we consider the expansions b k (x) = Yl\^i Cj^jfa) aR d b k /(x) = Yl\tli dj<fj(x), 
and observing that b k >(x) = (— l)^ +1 b k (x) and ^|A|-j( x ) = ( — l)^<Pj( x ) we have 

dj = {<PjM = X>i(*)(-l) W+1 &*(*) = - EviAI-i^ 6 *^) = " C I A I^- ( 2 - 3 ) 

X X 

Therefore 

H^A -2d)b k ,\\ 2 = Y \d j \ 2 (e j -2d) 2 = Y |c| A |_ j | 2 (2d-e| A |_ J -) 2 = \\(h A - 2d)b k f , (2.4) 

3 3 

and 

(bk',(h A - 2d)b k >) = Y \ d 3\ 2 ( e 3 ~ 2d ) = S ^Z l c |A|-il 2 ( 2d ~~ e |A|-j) = ~( b k, (h A - 2d)b k ). 

j 3 

(2.5) 

where the identity ej + e\A\-j = 4g? was used. 
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2.2. A bound for \\(h\ — 2d)bk\\ 2 - Now, it remains to show that the first term in the 
r.h.s of ( |2.2| ) can not vanish for all k in the fermi surface. First, let us consider the case 
d=2. Since 

[(h A -2d)b k ](x) = -Y J h(x + e) = -e- lk - x £ £ e - ik < e+e '\ (2.6) 

e e:x+e£dAx+e+e't£A 

the absolute value will be given by a sum of exponentials over some of the second nearest 
neighbors x + e + e' . The following diagrams illustrate the real part of the terms associated 
with each site, for particular values of k. Note that the configuration defines which terms 
will be in the sum. If we can prove, for suitable values of k, that all the terms have positive 
(or negative) real part, we conclude that they are not canceled out by each other, and a 
lower bound is obtained. 

o-l o+l 
+lo o o — 1 — lo o o — 1 

— lo o mx o o — 1 +lo o mx o o+l 

— lo o o+l — lo o o — 1 

o-l o+l 

=(§,§) fc=(0,±7T) 

Figure 1. Diagram of second nearest neighbors. 

For x G A, let [Q x ]ij = Qx,ij = e') : e || i,e' || j, x + e G dA,x + e + e' A}. If 

trQr ■■ 

2^ 2^ 

l^|[(/ iA -2 ( i)6 fc J(x)| 2 > |i J ^[(/ lA _ 2( i)6 fc J(x)| 2 > 1, (2.7) 

i=l i=l 

where the sum is take over fcj G {(±|,±|),(±|,=F^)}. Therefore we can conclude 

\\(h A -2d)b ki \\ 2 > #{xeA,trQ x ^0}, (2.8) 

for hi = (±5,+^) or (±^,+^). The same kind of argument makes ( |2.8[ ) valid for d = 3 
and fcj = (±§, ±f , ±§) or some vector obtained by inversion of coordinates. 
On the other hand, if trQ x = and Q x ^ 0, \{h A — 2<i)6^ (x) | 2 ^ 1 and 

|| (h A - 2d)b k \\ 2 > #{x G A : trQ x = and Q x ^ 0}, (2.9) 

for fc = (0, ±7r). For d = 3, the analogous result would be 

\\{h A - 2d)b k f ^ \#{x G A : trQ x = and Q x / 0}, (2.10) 

for fcj = (0, §,7r) or some vector obtained by permutation of the coordinates. 
If there are no isolated sites in A, 

#{x G A : Q x ^ 0} = a|0A| , a ^ (2.11) 

The reason that we need not consider the case where some of the sites in A are isolated 
lies in the fact that there is always a configuration obtained by joining this site to a 
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larger cluster, preserving the boundary. We only need to show that the energy of the new 
configuration (A') is lower than the original one. 

If we have a cluster and a disjoint site, the hamiltonian can be written as a direct sum 



v 2d / 

where h\ is the hamiltonian for the cluster. Consider now the perturbed hamiltonian 

where v{\) = (0, . . . , 0, —A, 0, . . . , 0), such that h A = h(0) and h A > = h{l). We know that 
the sum of first N eigenvalues {En) is a concave function of the perturbation A. Also, 
there is a unitary transformation that takes A — > — A, which implies that each eigenvalue is 
an even function of A. Combining these two results, we see that the sum of the eigenvalues 
is a decreasing function of A, and E-a,JV ^ E A t N . So, from this point on we can consider 
A as a single cluster. 

For simplicity, the remainder of the proof will be presented for d = 2, but the method 
is clearly general for arbitrary d. The only difference lies in the choice of the vectors k{. 
We are going to consider the two following possible cases: 

• Case I: #{x G A : trQ x / 0} ^ § |0A| 

\\{h A -2d)b ki \\ 2 > f|9A|forfc i = (±f,±f)ar(±f,Ti). 

• Case II: #{x G A : Q x + and trQ x = 0} ^ § |0A| 

\\{h A - 2d)b k .\\ 2 > f |9A| for ki = (0,±vr). 
For each of the two cases we have 

Now we need to know how rapidly can f{k) vary. 
2.3. A bound for |Vj/(fc)|. 

Lemma 2.1. For f{k) defined by (|2.2| ), the j-component of the gradient is bounded by 

< 10ad 3 . (2.15) 



3 \dk\ 



Proof: First, we should write 



JL\(h K -2d)b k {x)\ 2 = -i (ei + e 2 -e3-e 4 ) j e- 4fc ( ei+e2 - e3 - e4 ), (2.16) 

ei,B2,es,e4 

where the sum is taken over the edges such that x + e\, x + &% G <9A and x + e\ + e%, x + 
e 2 + e4 ^ A. We can bound the expression in parenthesis by 4, and the number of terms 
by (2(f) 4 . Also, the number of sites where \{h\ — 2d)bk(x)\ 2 does not vanish is limited by 
a|9A|. Therefore 

\^-\\{h A -2d)b k \\ 2 \ < a2 6 d 4 |aA|. (2.17) 

The same kind of argument leads to 

, d 



{b k ,{h A -2d)b k )\ ^ 6adt\dA\, (2.18) 
and combining the two results we conclude the proof of the lemma. 
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So, for 

..la 1 

|A; "^ 1 ^ w^d?'W¥ = 10 W (2,19) 

we have 

/(^) a 1 

|dA| ^ 2^ ^ 2^3- (2 ' 20) 

The lemma used to bound the gradient of f(k) is useful in determining a result for any 
dimension. However, if we focus on determining a better coefficient for d = 2, for instance, 
we should improve inequality ( |2.17| ). Instead of using the bound (ei + e2 — — e^)j ^ 4, 
we can sum over all possible vectors, using the real value of the expression in parenthesis. 
The same can be done for (2.18). The lower bound obtained for the j-component of the 



gradient of f(k) is 7a|<9A|. Therefore, it turns out that ( [2.20| ) is valid for the extended 
region 

We should make a remark concerning the fact that we don't know in principle which 
value of ki is the right one for Case I. But since e k is invariant under inversion of co- 
ordinates, the result will be the same, regarless of the choice between the neighborhood 
around h = (±§,±§) or k % = (±§,=F§) 

2.4. The lower bound. We are ready now to calculate the boundary term on the lower 
bound for the ground state energy. First, we should recall that 



E A , N - |A|e(l/2) > -±- d f(2d-e k ) £ ttU|(&*, <Pj)\ 2 dk. (2.22) 

^ ' J j-. ej >2d ^ ' 

But we proved that 

£ K^;)l 2 > f( k )\ dA \ > ^\ dA \ ( 2 - 23 ) 

j:ej>2d 

in the neighborhood of k{. We are ready now to state the preliminary result for d = 2. 
Proposition 2.2. Let region I be the neighborhood of ki = (tt/2,tt/2) defined by e k < 



£f(1/2) and (2.21), and region II be the neighborhood of ki = (0, ±7r) defined in a similar 



way. A lower bound for the ground state energy at n = 1/2 is given by 

E A , N - |A|e(l/2) > l^.-ljnunjf (2d- e k )dk = a x {l/2)\dK\. (2.24) 



where ai(l/2) > 10 



-17 



A similar result holds for d = 3. The regions I and II will be defined as the vicinities 
of the vectors ki presented in the last section. Also, we should include a factor of 1/3, to 
take into account (|2.1C| ). 

3. D = 2: A BETTER RESULT 

Considering a diagram like in figure 1 for the vector ki = (k,ir — k), we see that if 
Q x ^ 0, trQ x = and cos 2k > 0, 

\{h A - 2d)b ki (x)\ ^ cos2/c, (3.1) 
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whereas if trQ x ^ and cos 2k < 

^2\(h A -2d)b kt (x)\ > -4cos2&, (3.2) 

where the sum is taken over k{ G {(k,ir — k),(—k,ir — k),(k,—ir + k), (— k, — ir + k)}. 
Therefore we can extend the region of integration, as shown in the figure below. The 
shape of the internal boundary curve is defined by the cos 2 (2k) dependence. The new 
lower bound will be given by 

E A , N - |A|e(l/2) > . _L mm J (2d - e k ) cos 2 (2k x )dk = a x (1/2) |0A|. (3.3) 

Calculating the integral we get ai(l/2) > 1(T 13 , which proves our main result. 




Figure 2. Fermi surface e k = 4 and extended regions of integration. 



4. The result for n < 1/2 

For simplicity, we presented the detailed proof for n = 1/2. We should stress, however, 
that the method is quite general, and can be used to obtain the lower bound for the 
boundary term for any density n. Taking n = JV/|A|, we have an inequality which is 
equivalent to ( |2.2| ): 

,/ , m2 ^ ||(/ia - ejy)6 fc || 2 (b k ,(h A - e N )b k ) 
2^ l( ^' 6fc)l ^ 8d* + Ad ' (41) 

j:ej>e N 

Again, if (b k , (h A — ejv)£>fc) < 0, we take the application k' = k + (tt, it, . . . , ir), and we 
have e k t = 4d - e k , \\(h A - e N )b k \\ 2 = \\(h A - e\ A \_ N )b k > || 2 and (by, (h A - e\ A \- N )b k >) = 
~(b k , (h A - e N )b k ). 

Suppose we can prove \\(h A — eN)b k \\ 2 ^ a|<9A| for some k in the fermi surface. 
Then, either £ J:ej>ejv |(v?;A)| 2 > a/\dA\, or £* ej > e|A| _ w \(<Pj, V)| 2 > a'\dA\, which 
means that the boundary contribution can be calculated for density n or 1 — n. But due 
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to particle-hole symmetry, the boundary term should be the same for the two densities. 
Therefore, the problem reduces to proving ||(/ia — eN)bk\\ 2 ^ a|<9A|. 

This observation, when combined with the lower bound for the minimum of \\(h\ — 



£N)bk\\ over the fermi surface obtained in |FLU| is enough to determine the lower bound 



for the boundary term. To obtain a better and explicit coefficient, however, one should 
proceed like in the last sections, find vectors ki such that the function f{k) cannot vanish 
for all of them in arbitrary configurations, and take the integrals over the neighborhoods of 
these points in fc-space. The difference know is that we also need to perform the integration 
over neighborhoods of the vectors k[, situated in the image fermi surfaces of density 1 — n. 
Then we take the minimum over all of these integrals, to determine the lower bound. 
Therefore, the choice of the vectors ki might depend on the density, but apart from that, 
the method is quite general. 

5. Conclusions 

We showed here how to derive the lower bound for the boundary term of the ground 
state energy of Falicov-Kimball model. The existence of the boundary term is important 
since the system will try to minimize the boundary (to some extent) in order to minimize 
energy. Therefore, a segregated phase, where electrons and classical particles try to occupy 
distinct regions of the lattice, is obtained. When contrasted to the half-filling case, where 
crystalline long range order is observed, it might mean that the model has a first order 
phase transition when varying the chemical potentials. 

Our coefficient for intermediate densities is small when compared to the upper bound 



obtained in [FLU]. This means that our energy is not very sensitive with respect to 
the boundary size. However, the strength of the method is that it provides an explicit 
coefficient. Also, a\{n) can indeed be much smaller than ct2(n), since the upper bound is 
saturated by configurations with isolated sites, whereas the lower bound is not. 

The author is indebted to E. H. Lieb and D. Ueltschi for helpful comments and sugges- 
tions. 
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